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We examine the ground-state properties of the one-dimensional Hubbard model at quarter filling 
with Coulomb interactions between nearest- neighbors Vi and next-nearest neighbors ¥2- Using the 
density-matrix renormalization group and exact diagonalization methods, we obtain for U = IQt 
three different phases in the V1-V2 plane: 2A;f- and 4fcF-charge-density-wave (CDW) and a broad 
metallic phase in-between. Assuming that the metal is a Tomonaga-Luttinger liquid (TLL), we 
calculate the TLL parameter Kp. It is largest when Vi and V2 are frustrated, and Kp — 0.25 at the 
boundaries between the metallic phase and each of the two CDW phases. 
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PACS numbers: 71.10.Fd 71.30.+h 71.10.Hf 

The low-energy physics of one-dimensional correlated 
electron system can be described by the Tomonaga- 
Luttinger liquid theory^. In general, the Green func- 
tion and various correlation functions show a power-law 
behavior as a function of momentum k and frequency 
uj. The decay of the correlation functions is determined 
by the so-called TLL parameter Kp which depends on 
the strength of the interactions. Numerous experiments 
as well as theoretical studies have been performed on 
a variety of systems in order to determine Kp for one- 
dimensional metallic systems. 

Most features of a recent angle-resolved photoemis- 
sion spectroscopy (ARPES) experiment for the quasi one- 
dimensional organic conductor TTF-TCNQ^ can be ex- 
plained by the photoemission spectral function of the 
one-dimensional Hubbard model with only on-site inter- 
action^. However, the Hubbard model yields Kp > 0.5, in 
contrast to Kp ~ 0.18 as deduced from angle- integrated 
data for the density of states- which suggest a power-law 
behavior of the density of states near the Fermi level, 
p{uj) oc w", where a = {Kp + Kp^ -2)/4: with a ~ 1. An- 
other example of a small Kp is the quasi one-dimcnsional 
organic conductor (TMTSF)2X, where X^PFg, ASFg, 
or CIO4. In this system, the parameter Kp was esti- 
mated to be Kp ~ 0.23 from the power-law dependence 
dominating in the higher-energy part of the optical con- 
ductivity, (t{lu) ^ 10^°" ^p^^^ where a is the order of the 
commcnsurability [a — I for half filling; a — 2 for quarter 
filling)^-^. This value is consistent with ARPES measure- 
ments^. TLL behavior was also suggested for the strongly 
anisotropic transition- metal oxide PrBa2Cu408; the pa- 



rameter was estimated to be Kr, 



0.24 from both the 



optical conductivity and the ARPES study of Zn-doped 
PrBa2Cu40a2:3. 

Common features of all the materials mentioned above 
are that the filling of the conduction band is near one 
quarter, n ~ O.S^SiLiSiii. Consequently, they are close 
to a charge-density-wave (CDW) instabilitjiiSiiiiiii^. In 
addition, the parameter Kp is rather small, which is in- 
compatible with the result from a simple one-dimensional 
Hubbard model and shows that the long-range Coulomb 
interactions between carriers is relevant. Therefore, we 



examine the one-dimensional Hubbard model with inter- 
actions to nearest and next-nearest neighbors, 



(1) 



where c] ^ {ci^a) is the creation (annihilation) operator of 

an electron with spin ^(^ti i) at site i and ni ,j — c[ ^Ci^a 
is the number operator, t is the hopping integral be- 
tween neighboring sites, U is the strength of the Hubbard 
interaction, and the charge frustrating interactions Vi 
and V2 determine the nearest-neighbor and next-nearest- 
neighbor Coulomb repulsion. We restrict ourselves to the 
case of quarter filling (n = 1/2), and strong coupling, 

U = lot. 




FIG. 1: Phase diagram of the one-dimensional t-U-V\-V2 
model for U = IQt at quarter filling. The curves show the 
contours of constant TLL parameter Kp. Bold lines represent 
the boundary of the metal-insulator (CDW) transition and 
the dotted line corresponds to V2/t = Vi/{2t). 
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So far several theoretical studies have been made on 
this and similar modelsi^^iSiiSii^iSSiSi. In the ground 
state there exist two CDW phases with wavenumbers 
q = 2fcF and Akp. Between the two CDW phases 
there appears a wide region of a vanishing charge gap 
which results from the geometrical frustration of the long- 
range Coulomb interactions. If we assumed the inter- 
site Coulomb repulsions to be inversely proportional to 
the inter-site distance, we would find Vi = 2V2, so that 
neither of the two CDW instabilities dominates in the 
atomic limit, t = 0. Hence, one can easily imagine that 
the phase diagram contains a metallic region as soon as a 
finite t is introduced. However, little is known about the 
physical properties of this metallic state. It is the purpose 
of this paper to obtain an accurate ground-state phase 
diagram of the t-U-Vi-V2 model, including the TLL pa- 
rameter Kp whereby we make the natural assumption 
that the metallic phase is a Tomonaga-Luttinger liquid. 
Our phase diagram is shown in Fig. ^ 

We apply the density-matrix renormalization group 
(DMRG) method for the calculation of the ground-state 
energy. We study chains with up to 256 sites with open 
boundary conditions and keep up to m = 2000 density- 
matrix eigenstates so that the maximum truncation error 
is about 10~^. We use periodic boundary conditions for 
the calculation of charge excitations with the Lanczos 
exact diagonalization technique. 

In order to determine the metallic region in the phase 
diagram, we calculate the two-particle charge gap, 

A(2)(L) = Eo{N^ + l,N^^ + l,L) (2) 
+EoiN^ 'l,Ni~- 1, L) - 2EoiN^,Ni,L) , 

where Eo{N-^ , Ni, L) denotes the ground-state energy of 
a chain of length L with N-^ spin-up electrons and Ni 
spin-down electrons. When pairing is absent, in the ther- 
modynamic limit the two-particle charge gap becomes 
twice the single-particle charge gap, lim^^oo Ac {L) — 

2 limi^oo Ac(i). We find that Ac '{L) decreases mono- 
tonically with increasing L, so that we can extrapolate it 
to the thermodynamic limit systematically by performing 
a polynomial fitting in 1/L. Note that we need relatively 
large system sizes to obtain an accurate phase bound- 

(2) 

ary. The extrapolated results Ac /t at U/t — 10 and 

Vi/t = 4 are shown in Fig. |2] 

(2) 

It is evident that Ac is finite for both small V2/t < 

0.66 and large ¥2/1 > 3.43, i.e., the system is insulat- 
f 2) 

ing, and Ac /t vanishes in a wide range of V2/t, i.e., 
0.66 < V2/t < 3.43, within the accuracy of the extrapo- 
lation (error smaller than lO^^t). We find similar situa- 
tions for other values of Vi /t as well, which demonstrates 
that a stable metallic phase indeed exists between two 
insulating phases as shown in Fig. ^ This result is con- 
sistent with other studie9i2ii2ii2i22iSi. 

In the Tomonaga-Luttinger phase, the dominant cor- 
relations are determined by a single parameter Kp, which 
is given by22, 23^24 

Kp = {Tr/2)n^KV„ (3) 




FIG. 2: Charge gap Ar/C^t) (squares) and TLL parame- 
ter Kp (circles) in the one-dimensional t-U-V\-V2 model at 
quarter filling for U/t = 10 and Vi/t = 4. Lines are guides to 
the eye. 
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FIG. 3: Scaling of the ground-state energy differences from 
Eq. S'Hd the excitation energy from Eq. @ as a function 
of 1/L at U/t = 10, Vi/t = 1 and V2/t = l. The lines indicate 
the polynomial fits by Eq. ® to the data. 

where k and Vc are the compressibility and the charge 
velocity, respectively. To estimate Kp numerically, we 
calculate the compressibility k from the charge gap, 

A(2)(L) = 4/(n2L«(L)), (4) 

and the velocity of the charge excitations Vc from 

Eip{N^,N^,L) - Eo{N^,N^,L) - 2tiv,{L)/L , (5) 

where Eip{N^, Ni, L) is the lowest excitation energy with 
momentum k — 271 / L and total spin 5 = for finite 
system size L. As shown in Fig.|31 both of the quantities 
are monotonous functions of 1/L for all parameter sets 
and can be fitted to a polynomial function 

f{L)^aiL-^ +a2L-^ + ■■■ . (6) 

Therefore, we obtain 

K = lim k(L) = 16/a^'"^ , (7) 

= lim Veil) = a'^^''^/2Tr , (8) 

L — >-oo 
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so that the TLL parameter can be estimated as 



(9) 



in the thermodynamic hmit. For the compressibihty we 
can use DMRG for up to 128 sites. The extrapolation 
is very well behaved, and a reliable extrapolation could 
have been obtained from the results for much smaller 
systems. This makes us confident that the extrapola- 
tion for the charge velocity is meaningful despite the fact 
that exact diagonalization is limited to L < 20. A com- 
parison with the exact K'^^^'^'^ for the one-dimensional 
Hubbard model (Vi = V2 = 0)2^ shows that relative er- 



rors li^DMRG 



j>^exact 
P 



/^exact g^j.g ^jgiQ^ finite 
Vi and/or V2, our results are consistent with other esti- 
mates using the charge structure factori^^iSSiSi 

In Fig. 13 the parameter Kp is plotted as a function 
of ¥2/1 at fixed U/t = 10 and Vi/t = 4. We see that 
Kp > 0.25 when the charge gap is zero and Kp = 0.25 at 
the critical points. We have confirmed numerically that 
Kp is always 0.25 at the CDW critical points for all fi- 
nite values of U . Kp reaches its maximum value around 
Vi — 2V2, i.e., the density-density correlations decay 
most rapidly when Vi and V2 maximally frustrate each 
other. In general, long-range Coulomb repulsions are ex- 
pected to suppress the value of Kp. This is consistent 
with our results because Kp decreases when the values of 
Vi and V2 deviate from the relation Vi ~ 2V2, whereby 
the effective interaction strength increases. Apparently, 
the line Vi = 2V2 goes along the ridges of the contour 
line of Kp in Fig. This has been already suggested in 
the spinless fermion case and similar modelsiSiSi. 
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FIG. 4: Charge gap for spinless fermions at Vi/t = 20 (upper 
part) and Vi/t — 40 (lower part). The charge gap vanishes 
around V^/t = Vi/{2t) — 0.6 in both cases. The inset shows 
the charge gap for V2/t — 19.4 as a function of 1/L for system 
sizes up to L = 256. 



In the limit of large values of Vi and V2 {Vi,V2 < 
U), the boundary between the ordered phases and the 
metallic phase shrinks. In Fig. 01 we show the charge 
gap as a function of V2/t at Vi/t = 20 and Vi/t — 40 
for spinless fermions (U — 00). It demonstrates that a 



narrow but stable metallic regime exists between the two 

insulating phases. For large Vi and V2 it is very difficult 

(2) 

to extrapolate Ac reliably from exact diagonalization 
data (L < 20), and DMRG (L = 0(200)) must be used 
in the extrapolation. For example. Fig. 0] shows that 
Vi/t = 40 at V2/t = 20 is a CDW insulator whereas it 
has been assigned a 'non-Tonionaga-Luttinger metal' in 
Ref. I23 Even with DMRG it is very diflicult to decide 
whether or not the metallic region shrinks to zero at some 
tri-critical point, as proposed in Ref. |^ If it exists it 
is far beyond the values quoted previousljiSi. As seen 
from Fig. 0] there is a metallic region for Vi/t — 40 in 
the vicinity of V2/t = 19.4. In the surrounding of this 
point, the gap nicely scales to zero as a function of inverse 
system size, as seen from the inset to Fig. 0] Our results 
indicate that the metallic phase appears below the line 
V2 = yL/2, around V2 ^Vi/2- 0.6t. We speculate that 
there is a metallic phase between the two CDW phases 
for all finite Vi and V2, i.e., there is no tri-critical points 
in the V1-V2 phase diagram. 

In order to study the spin degrees of freedom, we cal- 
culate the spin gap for system size L, defined by 

A,{L) = Eo{N^ + l,N_l-l,L)-Eo{N^,N^,L). (10) 

As in the case of charge gap, the extrapolation to the 
thermodynamic limit is straightforward. In the 2kp- 
CDW phase, wc find that Ag is always finite because 
the system contains separated spin singlet pairs. For 
fixed Vi, As increases as a function of V2 and eventu- 
ally saturates at Ag = 4t'^/{U — Vi) in the limit V2 00. 
This is readily understood because, for U,V2 3> Vi, the 
system can be mapped to an effective spin Hamiltonian 

H = J2fcp E^ 84, • S4.+I with J2fcp = 4tV(i7 -Vi). This 

model is trivial and the spin gap is the energy differ- 
ence between the singlet and triplet state at each bond, 

A, = uyiu-Vi). 

In the 4fcF-CDW phase, we find that Ag is always zero, 
because a charged site and a vacant site come alternately 
and the spin degrees of freedom can be described in 
terms of a one-dimensional uniform Heisenberg model. In 
fact, for U > Vi/2 ^ V2, the effective spin Hamiltonian 
can be written as H = J^k^ X^i • 824+2 with J4fep — 

4 {e/{Vi - 2V2)f [ll{U - V2) + 2/{U - 2V2)]. This ef- 
fective Heisenberg model displays gapless spin excita- 
tions, in agreement with our numerical results. 

Lastly, in Fig. |5l we plot the TLL parameter Kp as 
a function of Vi/t for fixed Vi = 2V2 and several values 
of U. For large J7, Kp decreases as a function of Vi/t 
and eventually crosses Kp = 0.25 at some finite value of 
Vi.c. As shown above, Kp > 0.25 for a metallic phase, so 
that the Tomonaga-Luttinger liquid turns into the Akp- 
CDW at Vi^c- For small U, e.g., U = 2t in Fig. Kp 
decreases as a function of Vi , displays a minimum around 
Vi/t « 0{U/t) with -fCp, mill > 0.25, and increases again. 
This results from the fact that Vi overcomes the Hub- 
bard interaction U and electrons with opposite spin gain 
energy from on-site pairing. Eventually, Kp can become 
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larger than unity and superconducting correlations are 
dominant. As seen from Fig. \S\ the ground-state phase 
diagram in Fig. ^ is representative for all U ^ At when 
a superconducting phase does not interfere. The ridge of 
Kp goes along the line Vi — 2V2 and Kp decreases as the 
value of Vi and V2 deviate from this line. 




V^/t(=2VJt) 



FIG. 5: TLL parameter Kp of the t-U-V\-V2 model at quar- 
ter filling as a function of V\/t along a line Vi = 2V2 for 
U = 2t (circles), U = lOt (squares) and U = 00 (diamonds). 
The lines are guides to the eye. 

In conclusion, we obtained an accurate ground-state 



phase diagram of the t-U-Vi-V2 model at quarter filling 
using DMRG and exact diagonalization methods. For 
intermediate to large Hubbard interaction U ^ At, the 
system has CDW phases with q — 2kp and Akp be- 
tween which there appears a broad region of a Tomonaga- 
Luttinger liquid. Because of the geometrical frustration 
of the long-range Coulomb interactions, Kp is maximum 
around V2 = Vi/2. It is smallest at the phase boundaries, 
Kp = 0.25, which appears to be universal for transitions 
between the metallic and the CDW phasesSi. 
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